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Abstract. - We show how a qubit can be fault-tolerantly encoded in the infinite-dimensional 
Hilbert space of an optical mode. The scheme is efficient and realizable with present tech- 
nologies. In fact, it involves two travelling optical modes coupled by a cross-Kerr interaction, 
initially prepared in coherent states, one of which is much more intense than the other. At 
the exit of the Kerr medium, the weak mode is subject to a homodyne measurement and a 
quantum codeword is conditionally generated in the quantum fluctuations of the intense mode. 



Introduction. Quantum Information has recently entered the domain of continuous 
variable (CV) systems [1,2] and Quantum Error Correction theory in particular has been 
extended to the CV framework [3,4]. In such a framework, the main damage induced by 
decoherence is typically a small diffusion in the two canonical coordinates characterizing the 
CV system, which may represent the position and momentum of a particle or two orthogonal 
field quadratures of an optical mode. Gottesman et at have shown that these effects are best 
fought by encoding quantum information within appropriate shift-resistant quantum codes 
of the CV system [4]. The relevant aspect of the proposal of Ref. [4] is that, using qubits 
encoded in such a way, a universal set of fault-tolerant quantum gates can be implemented 
using only "cheap" standard devices which, in the case of a CV optical mode, correspond to 
beam splitters, phase shifters, squeezers, homodyne and photodetectors. In such a proposal, 
all the difficulties are confined to the preparation of the encoded states which, in the ideal 
limit, are superpositions of an infinite number of infinitely squeezed states and require the 
use of particular nonlinear interactions for their generation. The preparation of these shift- 
resistant quantum codes represents the main drawback of the Gottesman et al. proposal for 
universal fault-tolerant quantum computation in a CV system and therefore it is crucial to 
find an experimentally viable way to generate the quantum codewords employing optical CV 
modes. A first attempt has been based on a conditional scheme in which the CV system is 
coupled to a single qubit and which could be implemented on a single trapped ion [5]. Here, 
elaborating on a suggestion in [4] of using a nonlinear Hamiltonian of ponderomotive type, 
we present an all-optical scheme involving two modes interacting in a crossed-Kerr medium. 
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Ideal and approximate codewords. - Ref. [4] shows that a single qubit, whose generic state 
can be written as cq |0) + c\ |1), can be encoded in an oscillator with canonical coordinates 
q,p (with [q,p] = i) so that the corresponding logical state Co|0) + ci|l) is easily correctable 
with respect to small diffusion errors in both position and momentum. The codewords |0) 
and |1) are the simultaneous eigenstates with eigenvalue +1 of the two displacement operators 
T q (29) = exp {— 2ip9} and T p (2ir9~ 1 ) — exp {2iirq/9} (the so called stabilizer generators of the 
code [6]) with arbitrary 9 £ R. The two codewords arc a coherent superposition of infinitely 
squeezed states both in position and momentum, that is, 

+ 00 +OG 

]6)= Yl W = 29s)= Yl \p = *o- x s) (i) 

S — — 00 S— — OQ 

il)=f,(0)jO)= Y W = 29s + 9)= Y (-iy\p = nO-'s). (2) 

S— — OC s— — oo 

They have a comb- like structure with spacing 29 in the g-space and tt/9 in the p-space and 
are displaced with respect to each other by 9 along the q-axis (this is also responsible for the 
alternating signs in the p-representation, see Fig. ^| . 

In a CV system, errors and decoherence manifest themselves through random shifts of the 
two quadratures q and p. The present coding is designed for the case when these shifts are never 
too large, i.e., the error is correctable only if the shift is below a certain threshold. In fact, the 
recovery process is realized by measuring the stabilizer generators T q (29), which is equivalent 
to measure p (mod7r# -1 ) and therefore detects momentum shifts, and T p (2tt6 ), which is 
equivalent to measure q (mod 9) and detects position shifts. Then the error is corrected by 
displacing the state to the nearest peak, that is, to the nearest multiple of ir9~ 1 along the 
p-axis and to the nearest multiple of 9 along the g-axis. Therefore the recovery is successful 
only when the momentum shift Ap satisfies \Ap\ < ir9~ 1 /2 and the position shift Aq satisfies 
\Aq\ < 9/2, because otherwise the displacement would recover the wrong codeword, corrupting 
therefore the encoded information. An alternative choice for the codewords is the basis formed 
by the two superposition states |±) = (|0)±|l))/v2, which are used in [4] for the error recovery 
in the p-space because they have the same displaced comb-like structure of |0) and |1) in the 
p -space r ather than in the g-space (see Fig. 1). Practically, the two sets of codewords {|0), |1)}, 
{|+), |— )} are needed to measure the stabilizers T q (29) and T p (2ir9~ 1 ) respectively [4]. 

In real experiments one can only create finite squeezing and, therefore, only approximate 
codewords, |I) ~ jl), jo) = f q (fi)\T) and ]±) = (|0) ± \1))/Af± withjV± = {2(l±Re[(0 |Tj])} 1/2 . 
In this practical situation, even if we had only small position and momentum shifts satisfying 
the above conditions for a correct recovery, one would still have errors due to the use of 
approximate codewords possessing a nonzero overlap (0 |1) and (+ |— ). These intrinsic errors 
manifest themselves once the measurement of a stabilizer is performed. The quality of an 
experimental implementation of the error correction scheme of [4] is therefore established by 
an intrinsic error probability, which can be considered as the largest error probability among 
those associated with the four approximate codewords |0), |1), |±). 

The all- optical implementation scheme. - Our proposal for encoding a qubit in a CV 
optical mode is schematically described in Fig. [21 We consider two radiation modes a and 
b ([a,a^] = [b,b^] — 1) impinging on a nonlinear medium where they interact through the 
cross-Kcrr effect, described by the interaction Hamiltonian [7] 

H = HxtfaPb, (3) 
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Fig. 1 - Schematic representation of the position and momentum wave functions of the four ideal 
codewords |0), |1), |+) and |— ). 



where \ is the cross-Kerr susceptibility of the medium. The two modes are both prepared 
in coherent states, \a) and \f3), with (3 3> a (we take a and (3 real from now on). We are 
interested in the fluctuations of the b mode around the classical intense coherent state \f3), and 
therefore, redefining the b mode through the displacement b — > b + (3 and neglecting higher 
order contributions, the interaction Hamiltonian of Eq. can be rewritten as 

H = hx/3a?a(b + P) = -hktfa(b + P). (4) 



->n+AAM 

->[H 

Fig. 2 - Schematic description of the all-optical CV encoding scheme. 

If the losses within the Kerr medium are negligible the system dynamics is unitary and 
one has that the two-mode state after an interaction time t within the Kerr medium is given 

by 

\m = 7T- 1 / d 2 7 e-( Q2+ '^)/ 2 J2 ^r= { <»P{»**[<7* " nkt/2]} \n) a | 7 ) fc (5) 

where \n) a and |7) fc respectively represent a number state for a and a coherent state for b. 

Just at the exit of the Kerr medium, mode a is subject to a homodyne measurement of 
the field quadrature x = (a + a))/^/2 with measurement result x. As a consequence, the 
two modes are disentangled, with mode a in the corresponding eigenstate \x) a of x and mode 
b in the state conditioned to the measurement result x, which, introducing the quadratures 
q = (b + &T)/\/2 and p = (b — V<)/i^/2, can be described by the spatial and momentum wave 
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functions 

oc 

(p(q; t, a, x) — 7r~ 1 / 2 .A/~y^ (j, n (a,x) exp[— \q 2 + innrq] (6) 

n=0 

oo 

ip(p;r,a,x) = tt~ 1/2 N'^ /i„(a, x) exp[-|(p - 7rnr) 2 ], (7) 

where TV = tt 1 /^^,^ /w^e"* («-«') 2 r 2 } -i/ 2j = ^ x ) e -(" 2 +* 2 )/2 ; pii(a7x) = 

a n H n (x)/2 n / 2 nl, with H n [x) the Hermite polynomials, and we have introduced a dimension- 
less, scaled, interaction time r = \Z2kt/n. Our proposal consists in taking this conditional 
state as the approximate codeword |1), i.e., f(q] t, a, x) = (q\l) and ^(p;r,a,x) = (p 1 1 ) , with 
the following correspondence of parameters 9 «-> l/2r. The other approximate codeword state 
|0) is then obtained by applying a displacement, i.e., |0) = T g (l/2r)|l), while |±) are obtained 
using the fault-tolerant schemes discussed in [4]. 

Evaluating the performance of the encoding scheme . - From the wave functions of the 
encoded states shown in Fig. 1 it is easy to define and calculate the intrinsic error probabilities. 
In fact, an intrinsic error occurs when: i) given |1), the measurement of q (mod 9) gives 
a result inside one of the error regions R\ = [(2s — 1/2)0, (2s + 1/2)6'] for integer s; ii) 
given |0), the measurement of q (mod 9) gives a result inside one of the error regions R? s = 
[(2s + 1/2)0, (2s + 3/2)0]; Hi) given |— ), the measurement of p (mod n/9) gives a result inside 
one of the error regions Rj = [(2s— l/2)(7r/0), (2s+l/2)(7r/0)]; iv) given |+), the measurement 
of p (mod 7r/0) gives a result inside one of the error regions Rf = [(2s + l/2)(7r/0), (2s + 
3/2)(ir/9)]. Therefore, the corresponding intrinsic error probabilities are given by n 9 (r, a, x) = 

Yil=-oo i/fli I (q l 1 )! 2 ^, for both cases and ") and b y n ± = Ei^oo Jr± I (p l±)l 2d P for 

cases Hi) and iv). As discussed above, we can take as the overall intrinsic error probability 
of the proposed approximate encoding the largest error probability, that is, n max (r, a,x) = 
max[n g (r, a,x),n_(r, a,x),n + (r, a,x)]. It can be seen that this quantity is a decreasing 
function of the scaled interaction time r, at least for sufficiently large values of r. It is therefore 
convenient to choose a large r (which means using a large enough value of the coherent 
amplitude /3) and simply consider the asymptotic value Tl^ x (a,x) = Hm T ^ oo n max (r, a, x). 
The threshold value r t ^ of the scaled interaction time beyond which the asymptotic values of 
the error probabilities are essentially reached can be easily evaluated; r th depends upon a and 
for 0.5 < a < 2.5 it is r t ^ ~ 2 and therefore easily achievable with a moderate value of /3. 

The asymptotic error probabilities n^°(a,x) and II^ (a,x) can be evaluated analytically 
and are given by 

UQO( , 1 , 2 E^o(~l) fc (2fc+l)-Vnfa +2 ( 2fc+ i) 

n 9 (a,x) = - + - (8) 

>. | _ fc=o P2k+i / q s 

Ll±{a,X) - 2V oo 7-2 \ I 9 ) 

and their behavior as a function of the homodyne measurement result x is shown in Fig. 

It is evident from Fig. that the lowest intrinsic error is achieved for measurement results 
x as close as possible to x = 0, which is the measurement result for which II^ ax (a, x) is 
minimum. The probability densities in q and p for the four approximate codewords |0), |1), 



S. Pirandola, S. Mancini, D. Vitali and P. Tombesi: Constructing finite dimensional codeswith optical continuous v. 




-6 -4 -2 2 4 6 




-6 -4 -2 2 4 6 -6 -4 -2 2 4 6 

x x 



Fig. 3 - (a) Asymptotic error probabilities ilj°(a, x) (solid line), TV^(a,x) (dashed line), Ilj^(a, x) 
(dotted line) versus the homodyne measurement result x and for a — 1.5. The error probability 
II^„(a, i) is given by the largest value between these three functions and its behavior for three 
different values of a is shown (solid line) in (b) (a = 0.5), (c) (a = 1.5) and (d) (a = 3). The dashed 
line in (b), (c) and (d) corresponds to the asymptotic probability density, P°°(x,a), of getting the 
result x. 



and |±) corresponding to the particular measurement result x = are shown in Fig.^J where 
the parameter values a = r = 2 have been chosen. A first qualitative comparison with the 
ideal codewords of Fig. 1, shows that the present scheme is able to generate well approximated 
codewords. 

It is also evident from Fig. [3] that for values of x of the order or larger than a the error 
probability becomes unacceptably large. Therefore the coding scheme proceeds by fixing an 
acceptable interval centered around x = for the homodyne result, say —a/z < x < a/z, 
with z > 1, and accepting the resulting approximate codeword |1) only if x lies within 
this interval. In such a case one has two relevant quantities: i) the success probability of 
the coding scheme P(z,a) — f^™j* P°°(x,a)dx, where P ca (a,x) = n^ 1 ^ 2 J2^=o ^ s tnc 
probability density of getting the result x in the homodyne measurement (referred to the 
asymptotic limit of large r); ii) the corresponding mean intrinsic error probability H(z, a) = 
f-a/z P°°{ x > a )P{ z i a ) _1 nmax( a i x)dx. The success probability P(z, a) gives the fraction of 
times a satisfactory approximate codeword is generated. The mean intrinsic error probability 
11(2:, a) is the most relevant quantity since it measures the effective quality of the generated 
codewords. The behavior of these two probabilities is plotted as a function of z for different 
values of a in Fig. [S] It is easy to see that for z ~ 0, which means accepting all measure- 
ment results, 11(0, a) > 1/2, i.e. the encoding is too bad and cannot be used for quantum 
information processing. It is also evident that the best results are achieved for vary large z, 
i.e., when only values very close to x = are accepted, where however the success probability 
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Fig. 4 - Probability densities of the approximate codewords |0), |1), |±) for a — r = 2 (9 — 1/4). 
In (a) the probability densities | (q |0}| 2 (solid line) and | (q |1)| 2 (dashed line) are shifted by 8 = 1/4 
as expected, while in (b) the corresponding probability densities | (p |0}| 2 and | (p |1)| 2 (solid line) 
coincide as expected, with spikes at -k8~ 1 s — 4tvs. In (c) the probability densities | (q |+)| 2 (solid line) 
and | (q | — )| 2 (dashed line) differ only slightly, with spikes at 9s — s/4 while in (d) their momentum 
probability densities | {p |+)| 2 (solid line) and | (p \ — )| 2 (dashed line) are approximately shifted by 
TY0- 1 =4tt. 



becomes lower. In general the value of the coherent amplitude a has not to be either too large 
or too small. In fact, if it is too small, the intrinsic error probability is always too large (see 
Fig. [S] for a = 1), while if it is too large, U(z, a) stays around the value H(z,a) = 1/4 for a 
large interval of z (see Fig. [S]for a = 3,4, 5) due to the oscillating behavior of n^ ax (a, x) (see 
Figs. Etc) andOJd)). This fact implies that one has to use very large values of z to achieve 
an acceptable intrinsic error, corresponding however to a very low success probability. We 
see from Fig. that one has a good compromise between success probability and intrinsic 
error probability for intermediate a, such as, for example, 1.5 < a < 2. As an example, when 
a = 2 and z ~ 27 one has a mean intrinsic error probability 11(27, 2) ~ 1.0% and a success 
probability P(27,2) - 1.7%. 

Conclusions. - We have presented an all-optical conditional scheme able to generate the 
quantum codewords of the CV encoding scheme proposed in Ref . [4] for fault tolerant quantum 
computation. The generation of the encoded states is the most difficult part of that proposal 
and necessarily involves a nonlinear interaction. The present scheme involves two travelling 
beams initially prepared in coherent states interacting in a nonlinear cross-Kerr medium and 
is feasible with present technologies. 
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